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Abstract: In this article we focus on a semiclassical Schrodinger equation with matrix-valued potential 
presenting a symmetric conjoint crossing of three eigenvalues. The potential we consider is well-known in the 
chemical literature as a pseudo Jahn-Teller potential. We analyze the energy transfers which occur between 
the three modes in terms of Wigner measures. 



1 Introduction 

We consider the following Schrodinger equation 

1 ^-(g,0)=V'g(9)ei2(R2,c3), 

where e > is the semi-classical parameter and where the matrix- valued potential V{q) — Vpjriq) 
is 

/ qi (72/\/2 \ 

V<Z e R^ VpjTiq) - -qi (72/V2 . (2) 

V q2/V2 / 

Therefore, the potential V is subquadratic and the operator —^Aq + V{q) is essentially self-adjoint 
on the Schwartz space. Such systems naturally arise in quantum chemistry when studying the 
dynamics of large molecules in the frame of Born-Oppcnheimer approximation. The parameter e is 
small because it encounters for the ratio of the mass of the electron and of the averaged mass of the 
nuclei (see [T] or [TH] for the derivation of such equations). The solution ip^{t) does not have any 
direct physical interpretation but only quadratic functions of it: for example, for scalar equations, 
the position density \ip'^{q, i)P gives the probability of finding the nuclei in the configuration q E R^ 
at time t. Therefore, one is interested in the asymptotic description for the time evolution of the 
matrix- valued Wigner function of tp'^ (q, t) , 

W^^P' {t)){q, p) ^ i27r)-'' [ ij/{q~^v,t)^^"{q+^v,t) e'^-Pdv 

which plays the role of a generalized probability density on phase space. Indeed, the action of the 
Wigner function against compactly supported smooth test functions a G {'R'^ ,C^^^) is simply 
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expressed in terms of the semi-classical pseudodifferential operator with symbol a, which is defined 
forVei2(R2,C3) by 

op^ (a)V'(g) = {2Trey J a etP-(9-"V(«) dv dp. 

We then have tr / W'^{^l;){q,p)a{q,p)dqdp = (op^ (a)'0 7 i^)L'^CR? c^) ■ example, one recovers 
the position density by — tr VV^ {tp){q,p) dp. We are concerned here with the description 

JR2 

of the weak limits of the Wigner transform, these distributions are matrix-valued measure called 
Wigner measure or semi-classical measures (see [IS], [T7] or [H])- Besides, these matrix- valued 
measures fj, = (/iij) are positive in the sense that their diagonal elements fiij are positive Radon 
measures and their off-diagonal elements /Xi j for i 7^ j are absolutely continuous with respect to 
and fj.jj. 

The potential Vpjriq) arises in models of quantum chemistry: it is the simplest pseudo Jahn- 
Teller Hamiltonian as introduced in 8J Chapter 10 (see also [30] or [10] for example). Its main feature 



is that its eigenvalues are y qf + q2J ^y^i +92 ^^'^ ^- They are symmetric: the eigenvalue is 
exactly the half sum of the two other eigenvalues. Besides, they are of multiplicity 1 as long as 
{ill (12) 7^ (OjO) E^nd they simultaneously cross on the point {qi = q2 = 0}. We will say that V 
presents a conjoint symmetric crossing in q = (0, 0) = 0. The appellation of pseudo Jahn- Teller 
potential is a reference to what is called Jahn- Teller potentials in quantum chemistry, namely the 2 

by 2 potential of the form ( '^^ 92 j -v^hich numerous mathematical works have been devoted. 

V I2 -qi J 

Then, another simplest realization of a symmetric conjoint crossing is the matrix 

(91 92 
92 -91 


Both matrices Vjt and VpjT (defined in ([2])) have the same symmetric eigenvalues and the behavior 
of solutions to ^ with the Jahn- Teller potential Vjt is well understood. Indeed, in that situation, 
the third mode evolves at leading order independently from the two other ones which interacts 
following the well-known process of conical intersections (see [22], [3T], [12], [lU], [13] for example). 
However, the situation is different for the pseudo Jahn- Teller potential Vpjriq), we will see that the 
three modes will interact altogether; note also that for this potential, the three spectral projectors 
are singular on {gi = q2 = 0} while Vjt has a smooth eigenprojector (for the mode 0). While 
an important literature have been devoted to two by two crossings, there are very few results on 
conjoint crossing of three eigenvalues. The aim of this work is to analyze on a simple model the 
mechanism of a symmetric crossing of three eigenvalues. 

We shall describe how the three modes interact together for the pseudo Jahn- Teller potential 
Vpjriq) by studying the branching of Wigner measures on the crossing set {qi — q2 — 0}- We explain 
our result in the following section. The main ingredient of the proof are a normal form result given 
in Section [3] and a scattering theorem in Section UJ The normal form reduces the analysis of the 
crossing to the study of a system very close to the system of three ordinary differential equations 

/ s z/V2 \ 

iedsu+ -s z/V2 = 0. (3) 
V z/V2 z/V2 / 
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The behavior of solutions of this system as e goes to is analyzed in Section [H This system has 
been already studied by physicists ([5] and [31); with their arguments, one can find the same result 
than ours, however, we give here a different and slightly shorter proof. These two results allow to 
understand the evolution of Wigner measures of families of solutions to ^ in Section [51 

This paper is the first step in the understanding of symmetric conjoint crossings. For applications 
in quantum chemistry, it would be interesting to be able to treat general potential V{q) presenting 
a conjoint crossing of symmetric eigenvalues on a codimension 2 submanifold. Then, it is likely 
that the potentials Vjt and Vpjt will play a crucial role as toys model. This motivates the work 
performed in this paper. Besides, the application to quantum chemistry that we have in mind 
concerns the elaboration of algorithms modelizing the evolution of the Wigner transform of families 
of solutions of Schrodinger equation ([1]). These algorithms have been extensively developed since 
the 70's and the pioneer work of TuUy and Preston (see [29]). The algorithms of [24] and [15] which 
are the first one presenting a mathematical proof relies on normal forms results of [Ij, [S] and [T^] 
and on the precise analysis of the propagation of Wigner measures through the crossing performed 
in [14] and [11] . Therefore, the result of this article concerning the branching of Wigner measures 
has important consequences for numerics in quantum chemistry: one can construct an algorithm 
for the potential Vpjx {q) exactly as in [24] and its convergence is guaranteed by the Theorem [1] in 
the same manner than the algorithm of [H] relies on the branching theorem of [13]. The diagrams 
illustrating the statement of the main result in the next section are obtained that way. 



2 Main result 

Let us begin with some notations. For £ e {0,+l,— 1}, we denote by n^(g) the eigenprojector 
associated with £ \ q \ and we set 

X\q,p,r)=T+^-^+£\q\. 

The indices £ € {0, +1, —1} will be sometimes shortened into ± whenever £ — ±1. We consider the 
Hamiltonian vector fields 

H^. =dr+p-V,~£^-Vp, £e{0,+l,-l} 
and their integral curves pi which are called classical trajectories. They are of the form 
pi^ {ql,s + to,pi,To) with ql^pi, pi^^-rY\ ^'^^ Po = (9o, to,Po, tq). 

ml 

Because of the singularity in q — 0, something has to be said for the ± mode: it is proved in |13j 
that under the assumption 

Po ^ 0, (4) 

there exists a unique curve {qf,pf) solution of the system of ordinary equations above and such 
that q^ — qo and — pq. Moreover, if — 0, the + trajectory (resp. — trajectory) issued 
from {qo,Pa) smoothly continues the — trajectory (resp. + trajectory) arriving at {qo,Po)- Finally, 
we emphasize a specific feature of these trajectories: there exists an hypersurface containing all the 
curves which pass through the crossing set, namely the singular hypersurface 

I = {qAp^O} 



3 



where for two vectors z = (zi, z-i) and z' = (zj^, z'^j of R^, z l\ z' denotes their wedge product: 

Z l\ z' — Z\z'^ — Z2z'i. 

Observe that / is a smooth hypersurface close to points po = (90 : io,Po, To) such that po 7^ 0. Note 
also that this assumption implies that the classical trajectories are transverse to the crossing set. It 
is under this assumption that most results on crossings are obtained (for example in [19j . |20j . [13] 
or [Uj; for a result in case of tangency see j9]). 

We work in space-time variables and call crossing set the points of the phase space where there is 
an eigenvalue crossing namely the set 

^=|r + ^ = 0, g = o| ={^? = 0}nE, 

where S denotes the characteristic set: 

S = E+ U U 11° with = {\\q,p, r) = 0} 

being the energy surface for the mode The Wigner measures of families (■0^) are supported on 
the energy surfaces. More precisely, it is proved in [18 that outside the crossing set, one has 

/i = /i+ n+ + + ^" n° 

where for i E {0, +1,— 1}, the measures fi^ are scalar positive Radon measure such that for all 
aGC„-(R^C), 

(op^ (a^'^)V'^ V'')i2(R,3^c3) ^ "('?'^'-P''^)d/(9,^:^',T)■ 
Besides, the measures /x^ are supported on and invariant through the Hamiltonian flow associated 
with Hxi for aU £ € {0,-1, +1}. 

It is well known that the transitions between two modes cannot be described only in terms of Wigner 
measures (see |13j). The phenomenon of the transitions is more intricate and requires a second level 
of observation. Indeed, the transitions are determined by the way the families ('ijj'^{t)) concentrate 
on the trajectories arriving at the crossing with respect to the scale y/e. For this reason, one uses 
two-scale Wigner measures associated with the hypersurface / because all the trajectories passing 
through S are contained in /. 

Following [2^ and [13], we consider observables a — a{q,t,p,T,ri) G C°°(R7 ,C^''^) which are com- 
pactly supported outside {p — 0} in the variables {q,t,p,T), uniformly with respect to the variable 
r] and coincide for I77I > Rq for some i?o > with a function homogeneous of degree in the variable 
rj. We will denote by A the vector space of such functions. Then, when given the hypersurface /, 
one associates with a the pseudodifferential operator 



op, 



/ (a) =op^ ^q,t,p,T, ^^^^ 



The change of variables 



iq,p) ^ {^/eq,^/ep) 

and Calderon-Vaillancourt Theorem show that this family of operator is uniformly bounded with 



respect to e on Lf^^CR"^). A two-scale Wigner measure associated with the concentration of the 
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family (■0^) on / is a positive Radon measure on the compactified normal bundle to / such that 



(op/ (a)V'^V'^)i2(R3,C3) 



a{q,t,p,T,'i])di^{q,t,p,T,ri) + I a ( q,t,p,T, ^ oo ) dfi{q,t,p,T). 



yj. 1 I J. I 11/ \ J- ' III I I / I \ J. ' 111 I \ . 

N(i) J I- \ \q^p 

We recall that the fibres of the normal bundle to / is obtained above a point p of / by quotienting 
the vector space Tp(r*R'^) by the tangent to /, Tpl. The fibre above / is a dimension 1 vector 
space Np(I) and one gets the fibre of the compactified bundle to / by adding two points at infinity; 
then, this fiber is isomorphic to R. The function a being homogeneous in the variable rj, it defines 
a function on Np{I). We point out that fj, above / is the projection of ly on /: 



tJ-{q,t,p,T)li ^ _iy{q,t,p,T,dri). 
Jr 



Since /i is determined outside / by the Wigner measures of the data (because of the invariance of 
/i^ through Hxt outside S), we focus on what happens above / and thus on the measure v. 
Similarly to Wigner measures, the two-scale Wigner measures satisfy localization and propagation 
properties: it is proved in (13j that outside S, 

1/ = n° + ty+ n+ + 1^- u- 

where for any £ S {0, +1, —1}, the scalar positive Radon measures i'^ are supported on and prop- 
agate along the linearized Hamiltonian fiow induced on N{I) by the Hamiltonian vector fields H^^t. 
We are now able to describe the branching of two scale Wigner measures close to a point po = 
(go, io,Po, To) such that po ^ 0. Indeed, since po 7^ 0, the classical trajectories are transverse to the 
crossing set S and for any I S {0, -1-1, —1}, the measures have traces on S. Let us denote by 
(resp. I'oui) the traces of on the in-going (resp. out-going) side of S. For expressing the link 
between the traces, we need to choose coordinates on N{I): for p = {q,t,p,T) £ I, we characterize 
the class in Np{I) of the vector 6p = {6q, St, Sp, St) of Tpl by the coordinate 

71 = Sq A p. 

These traces are linked as stated in the following theorem. 

Theorem 1 Above po = (qo, ^o,Po, 'Ha) G S with po ^ and under the condition that u^^, uf^ and 
are two by two mutually singular on {\rj\ < +oo}, one has 

(1-T)2 r2 2r(l-T) 

r2 (l-r)2 2T{1-T) I ( ^,r„ I (5) 

2r(l-T) 2r(l-T) (1-2T)2 







^out 





whe 



T(p,ry)=exp(-|J). (6) 



The coefficient T{p,r]) is very close to the Landau-Zener coefficient exp which appears for 

crossings of two modes (see [12], [SI], EO] and [13]). 

We could also have expressed the transitions through the crossing set by considering the two-scale 
Wigner measures associated with the concentration of {ip'^) on the sets Jf„ (resp. Jout) consisting of 
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all the trajectories for the mode I arriving to (resp. going out of) the crossing set for £ g {0, +1, —1}. 
These measures have traces on S that can be identified to the measures , .. The link is then 

in/ out 

given by ([5]). This point of view (which is the one of T3J) is more intrinsic and applies in situations 
where there is no trivial sets containing all the trajectories passing through the crossing. This is a 
real issue for general symmetric conjoint crossing (see [?]). 

We have reduced the dimension space to 2 but we could have assumed the dimension to be greater 
than 2 with a c?-dimensionnal Laplacian and a potential Vp,7t('?i, 92) depending only on qi and (72- 
Our result extends to this situation replacing / by {(^1,^2) A (91,(72) = 0} and p by (pi,P2) in the 
non-degeneracy condition @ and in the transition coefficient (O. It also extends to the situation 
where the Laplacian is replaced by a Fourier multiplier A{D). One then need to turn p into V A{jp) 
in (|4]) and ([6]). In that situation, the existence of a submanifold satisfying similar properties than / 
is non trivial. 

The proof of this theorem relies on a normal form result which is the subject of Section [3l One 
reduces the analysis of the crossing to a scattering problem which is studied in Section ID Provided 
these two results, the proof of Theorem [T] is performed in the last section following the method 
introduced in [131 (slightly modified to take into account the difficulties induced by the presence of 
the third mode). 

Before passing to the proof of the Theorem, we illustrate our result by some numerics. We choose 
a Gaussian initial data supported on the + level of the form 

= (e7r)-l/2exp(-2^|5 - ql\^ + ^ Po ■ (q - q'o)) e+{q) 

where 

e+(q) e Rann+(g), e - 10"^, po = (-l,0), = (^5, . 

The time evolves between t — and t — j and we calculate an approximate value of the population 
on each mode: 

n'{t)= [ \n'{qW{q,t)\^dq, ^e{-l,0,+l}. 
J-R? 

Following [24j and |15j . we perform a surface hopping algorithm: 

• we sample the Wigner function, 

• we propagate weighted points {qf,pf,w) along the classical trajectory until the distance to 
the gap {q — 0} is minimal, 

• at that point (gji , p J. , w) one generates new trajectories for each mode transporting part of 
the weight w: 

{qt,pt,a-T*)'w), {q-,p:,T*^w), {qlpl2T*{l-T*)w), s > s* 
where the result of Theorem [T] motivates the choice of the transition coefficient 



T* = exp 



2e\p[ 



,0 13 



Arguing as in one can prove the convergence of the algorithm. No estimation on the rate is 
available for the moment since it requires an improvement of the normal form result which is out of 
reach for the moment. In order to illustrate numerically the validity of the algorithm, we calculate 
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a reference solution based on a grid discretization and a Strang splitting scheme. This scheme is 
known to have a very fast rate of convergence (see [53], [T5] and [S]). Then, we compare the two 
outputs in the figure below: the initial population which was on the + level splits after passing once 
through the crossing. 




Figure 1: Evolution of the population on each level for surface hopping algorithm (full line) and for 
Strang splitting scheme (dotted line). 



3 Normal form 

In this section, we first prove a microlocal normal form which holds in the phase space close to a 
point of the crossing set. Then, in the second subsection, we write a two-microlocal normal form 
which is only valid in neighborhoods of size of the hypersurface /. 

3.1 The microlocal normal form 

The phase space has the structure of the cotangent space T*R'^. It is endowed with a symplectic 
structure given by the 2-form 

uj — dj = dp A dq + dr A dt 

where 7 is the Liouville 1-form 7 = pdq + rdt. A local change of sympletic coordinates k is called 
canonical transform and one can associate with it a unitary bounded operator of i^(R'^) compatible 
with the change of coordinates. Such an operator K is called Fourier integral operator associated 
with K. It satisfies Egorov's Theorem 

Va e C^CR^), K* op, (a) K =op, (a o k) + 0{e^) in C{L^(R^)). (7) 

We point out that the remainder is better than expected because we have chosen to use Weyl 
quantization. The reader will find in [57] a complete analysis of Fourier integral operators. 
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We consider ft a neighborhood of po = (9o , : Po 7 tq ) G S with po 7^ and we suppose fl small enough 
so that p ^ in n. Then, for £ G {0, +1,— 1}, we denote by j/„ (resp. Jout) sets consisting 
on all the classical trajectories for the mode i entering in (resp. ongoing of) S nV,. We point out 
that by the property of continuation of the trajectories, the sets (^J^ U J^^.^) n fl are submanifolds 
of T*R^ just as (J"„ U J°„f) n does. Our result is the following 

Theorem 2 Let po — (qo , to , Po , tq) G S be such that pa 7^ 0. Then, there exists a local canonical 
transform k from a neighborhood of po into some neighborhood of 0, 

K ■■ {q,t,p,T) {s,z,cr,C), k(po) = 0, 

a Fourier integral operator K associated with k and an invertible matrix-valued function A such that 
if Tp^ is a family of solutions to ([T]) for some initial data ipQ uniformly bounded in L^(R^), then 

= K* op, (A) V'" 

satisfies for all (p G C^{fl), 

op, (0) op, {-a + VpjT{s,a{C,a)zi)'^v' ^0{e) (8) 

in L'^{'R?) where a is smooth and does not vanish in fl. 
Moreover, 

k{I) = {zi = 0}, (9) 
k(J±) = {(7±s = 0, zi = 0, s<0}, (10) 
^(^o'^ut) = {^Ts = 0, zi = 0, s>0}. (11) 

Proof: The proof of this Theorem proceeds in two steps. First, we find the canonical transform k. 
Then, we use the properties of the matrix VpjT to construct A such that 



^ ( ( ^ + ^ ) Id + VpjT{qi,q2) ) A* 



2 

Then, we obtain ^ by the Egorov theorem ([7|) 



OK, = —(T Id + Vpjt(s, z) 



First step: the canonical transform. This first step relies on the analysis of the geometry of 
the crossing. We crucially use Proposition 6, p. 148 in [T3] . 

Proposition 3 There exist a local canonical transform 

K : {q,t,P,T) 1-^ {s,z,a,C) 

and non-vanishing smooth functions A and p such that 



o- A(p,t) 



m: 

2 



. = A(p,r)^.g, (12) 
zi = A(p, t)^(p,t) ill A q. 

Besides, one can choose ^ — C{p, t) such that [p, t) i-^ (C, a) is a diffeomorphism and 

)\s 



A(p,r)|5 = br^/^ (13) 
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Let us shortly describe the proof of Proposition 6 in [T3] . Observing that 

we obtain the existence of A such that s and cr satisfy the bracket condition {cr, s} ~ 1. Therefore 
we necessarily have (|13l) . Besides, for any function fi ~ t), one has 

Therefore, it is enough to find [i such that 

A(p, r)— • g , A(p, r)/i(p, r)— A g| = 0. 
Once is built, one completes the symplectic system (ct, s, zi) by Darboux Theorem. 

By the definition of S and /, one has niS) = {cr = s = zi = 0} and — {zi = 0}. Moreover, 
since A ^ (since p ^ in il), we have = {cr((T^ — — Zj) = 0}. In order to find the precise 
equations of k( Jj^) and k(J^j), we observe that 

<Jtlout)^{'^''-^\ ^1 = 0} (14) 
and we consider the Hamiltonian vector fields ll\± . We have outside S 

m 

Observing that if the trajectory {qf,pf) passes through at s = we have 

qf = spo + o{s) 

we deduce 

77;,±(p±) = 5t +p • V, T sgn(s)^ • V + o(I). 

bo I 

By the definition of a and s, we have 

Ha±s = A 

Therefore, we obtain 



\P\ 



on S. 



\Po\ 



This equation impHes that k{J^) C {cr + s = 0}, whence the equation of k{J^) by and 
dimension considerations. One argues similarly for Jg^^, and J^^. 

Second step: the gauge transform. We now transform the symbol (t + -l^j Id + V{q) in order 
to see the coordinates (s, z, cr, C) in the symbol. We first use the matrix 

M = 
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and observe that M = M* = M''^ and MV{q)M = W{q) with 

(0 qi (72 
91 
<Z2 

We then use the foUowing lemma which comes from a straightforward computation. 



Lemma 4 Consider the matrix R(p) = -p-r 

\P\ 

and we have 



We can now conclude the proof of Theorem [21 Set 



\p\ 





- 


-Pl 





P2 


-w 


(p_ 


\\p\ 



-Pl 



P{q,P,T) = ( T + ^ ) Id + V{q) and R{p) ^ MR{p)M. 



We have 

R{p)P{q,p,T)R{pr = (^r+M!^Id + l-(^^.g,^A 
Therefore, using the canonical transform of Proposition [31 we have 

X{p,T)R{p)P{q,p,T)R{py o K"^(s, z,a,Q = crid - V{s,a{p,T)zi). 
Let us consider now a Fourier integral operator U associated with k and set 

A{p,T)^X{p,Ty/^Rip). 

We have 

U op, [A] op, (P) op, {A*)U* =op, (aid - O^i)) + 0(e) 

in£(L2(R3)). □ 

3.2 The 2-microlocal normal form 

In this section, we want to ameliorate the normal form in open sets localized at a distance of order 
^/e of the hypersurface /. For simplicity, we set 

/10 0\ ,/00l\ 
J= 0-10 ,L=-= 1. (15) 

\oooy \ 1 1 J 

We shall take advantage of the fact that if a e C^(R^), 

op,-'^ (a) op, {-aid + Vpjt{s, a{aX)zi)) ^opj (a) op,-'^ (-aid + Vpjt(s, a(cr, C)\/e?7)) • 
We set 

Q = -ald + sJ + ^/ea{aX)vL, (16) 
Qo = -<jU + sJ + V^a{0,C)vL. (17) 

We prove the following result. 
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Proposition 5 For every ball B of R,,, there exist smooth matrix-valued functions C and C such 
that for all a e C§°(R*^ x B), 



op/ (a) op/ (Id + V^C) op, (Q)- op, (go) op/ (Id + V^C) 



cm 



0{e). 



Proof: We want to realize 

op/ (Id + ^/IC) op, (Q) =op, (Qo) op/ (Id + x/iC) + 0{e). 

For this, we successively consider the terms of order 1 and ^/i in the development of the preceeding 
equality by symbolic calculus. The terms of order are equal and the equality for the terms of order 
■v/i reads 

a(o-, C)r? L - a(0, C)?? L + a{C - C) + s{CJ - JC) = 0. 
We choose C = - '^C with 

/ 1 
C = f\ 1 
\ ^ 

so that the second equation writes in view of CJ = and JC = — *(CJ) = 

L = 0. 



Hence 



7?(a(a,C)-a(0,C))- V2a/ 

T] Q;(a,C) - q;(0,C) 



f{a,z,r]) = 



^/2 



This closes the proof of the Proposition. 



□ 



4 Scattering matrix for the pseudo Jahn- Teller model 

In this section, we analyze the differential system 



id.. 



s z/V^ \ 

-s z/V2 
z/V2 z/V2 / 



V- = Vpjt{s,z) V- 



where 2; is a non-zero complex parameter. We prove the following scattering result. 
Theorem 6 There exists (a~,a+,a'^) and (w~ , u)~^ , oj^) such that 

\zr 



vt{s) = exp 
ve{s) = exp 



'i{- + ^Ms\ 



a +0 { I as s ^ —00, 

uj*^ + { ] as s —^ +00. 



Besides 











= S{z) 











(18) 
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with 

\ -\/2e-*^/4f^5,(^) ^/2e'''/'^n{z)e{z) 2e"'^l^l'/2 - 1 



5(z) = 
and 



Note that the indices +, — and are not directly connected with the modes A^. However their 
connection will be clarified in Section! 



4.1 Providing two independent solutions with contour integrals 

Classical theorems give that the solutions form a tridimensional linear space of analytic functions. 
The parity block decomposition of Vpjt{st z) induces that the space of solutions split into two 
subspaces: a bidimensional one with v± even and ug odd, and a unidimensional one with v± odd 
and Vq even. For an even (respectively odd) analytic function w, there exists an analytic function 
/ such that v{s) = f{s^ /2) (respectively v{s) = sf{s^ /2)). Hence the system, in the case of the 
bidimensional subspace, 

dt Q 

a* I I /- I =Vpjt{1,z) I /- I . (19) 



1 + 2td, 




To solve ([T^. we will use the method of contour integral (see for instance ^ Chap. XV): we look 
for solution of the form 

f+{t) \ r ( F+{r) \ 

Ut) = / e"- F^{r) dr (20) 
fo{t) J \ Fo{r) J 

where A is an oriented path of C and the F^^s are fimctions holomorphic in a neighbourhood of A. 
Injecting ^ into (HH) we get 





and after integrating by parts 
















1] 














A 







I dT 




-T+ i 

-r-1 ^ I I F^(r) 1 dr. 
73 



Hence, it is enough to choose the i^*'s such that 




-T+1 





V2 \ 







-T - 1 


- j 


F-ir) 


z 

V2 


z 

V2 




\ Fo{r) 
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and A such that 



0. 



(21) 



The system is exactly solvable and gives, up to a global constant, 

i^[ln(r + l)-ln(r-l)] p ( \ 



V2- 



1 



z Fo{t) 
V2- 



1 



Now F^,{t) has three singularities, namely ±1, and F^,{t) is multivalued. More explicitly, consid- 
ering a closed contour, the value of F^{t) is changed according to the following rule: if the contour 
goes n± times anticlockwise around ±1 then ln(T=F 1) is changed into ln(r =F 1) +«27rn± and if it goes 
uq times anticlockwise around then ^/T is changed into (— 1)"°-\/t. If we choose a closed contour 
A that avoids the singularities of and solves (|2T|l . the corresponding function ([20|) is then well 
defined and is a solution of ([TO]) , the only thing that one has to check is that it does provide a non 
trivial solution. We will focus on the following contours: they enclose only the singularities ±1 and 
go once around each of them. There are exacly two of them depending on whether we stay below 
or above the singularity (in the first case the contour has the shape of a U and the second of a 
n). The idea is now to use residue arguments to get explicit formulas for the contribution of the 
different singularities. Multiplying by —y/2/\z\, we define our two solutions to be 




z 1 

V if^ 



dr 



(22) 



where the A*-^ and A*^ contours are decomposed according to Figure [2] and the determination of the 
two logarithms and the square root is such that it coincides with the usual definition (for positive 
real numbers) at the point M. 





Figure 2: decomposition of the A*-^ and A*^ contours 



For the sake of clarity, the points A, B, C and D (respectively A' and B') have been split but they 
have to be thought of as the same geometric point. To study separately the different pieces, we will 
first show the following residue results. 

Lemma 7 Let [3 he a non zero real number, w{t) an holomorphic function of t in a convex open 
neighborhood Q of a point a of the real axis and A an oriented path of Q with both ends in the upper 
half plane. Then, in the limit \s\ — > +oo 
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1. if A does not surround a, 



^ ^ J A \Jt - a 



dT 



= O (e-"^^ 



where S is the minimum of the imaginary part of the two endpoints 



Figure 3: non-surrounding and surrounding contours 



2. if A does surround a once and anticlockwise, 



g*V^e*'''"(^-")w(r)dr = -Aie'^^'^^e''^^^ sinh(7r/3)r(l+i/3)w(a)- 



1 + 



''w{t)- 



dr 



2%/2^e"/^u;(a)- 



1+0(4 



-a-2(3\n\s\ | p ^ ^ 

(23) 

. - . ^ , , (24) 

/A Vt- a 

where the determination of the logarithm and of the square root is such that it coincides with 
the usual definition for positive real numbers. Furthermore (|23p is locally uniform in [3. 

Proof: 

1. It follows from a crude estimate of the term e*""^ for the straight line path. 

2. We deform the contour into the one presented in Figure [3] where the two vertical lines have 
been split for a better understanding i.e. B and E (respectively C and D) are in fact the 
same geometric point. The radius of the small circle is chosen to be A(s)/s^ where A(s) — > 
as \s\ — *■ +0O. 

The pieces AB and EF are exponentially small by the first point. The piece BC is —e^^^ 
times (respectively equal to) the piece DE for the logarithm (respectively the square root). 
Let us compute the piece CD parameterizing the circle by a + ^^e*^ for — 37r/2 < 9 < n/2: 



2 ^g-2i/31n|s|gi/31nA(s) 



tt/2 



3ir/2 



■ Ms) ie 



= O 
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B 



V 



A 



D 




Figure 4: deformation of the surrounding contour 



/' 

Jc 



dT 



\/t — a 



e 2 



■37r/2 



o 



Finally, we compute the contribution of the piece DE parameterizing the vertical segment by 
a + 2^u for A(s)/2 <u< 



f 

Jd 



w{T)dT 



gi^ag-2il3 lii|s| 



gi/3i„2g-./3/22^^' e—u'f^w(a + 2^u]du 



-2i llu 



,ifiln2g-7r/3/22j 



gj-%-ag-2i|81n \s\ 



l\(s)/2 

){a) J e-'^u'f^du + O [\{s) + - 



gi/31n2g-./3/22^p(l + i(3)w{a) + O ( ^ + JL j , 



vt - a 



.s2<5/2 



\(s)/2 



e "u ^/^u; ( a + ] du 



S\ 
,1^5- a 



(a) e-^u-'/^du + O i^^/Ml) + ^ 



e 2 



-V2e^"/4M;(a)r ( - ) + O 



The result then follows from the choice A(s) = (respectively A(s) = 1/ s^) for the logarithm 
(respectively the square root) and the fact that r(l/2) = ^/^T. 



□ 



Let us now give explicitly the contributions in the case of the A*^ contour proceeding with 
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integration by parts and several applications of formula 



z 

2i—e 


-^|2|'/8 


\z\ 





















^l^p^g Sinh (^k|V4) / 
|z|V4 

(^klV4) 
|z|V4 



■\n\s\ 



0(^1 (25) 

(26) 
(27) 

The asymptotics for the A*-^ contour can be mostly deduce from the following considerations (up 
to exponentially small errors): 

• the contribution of the piece CD is the same as B' A' whereas the contribution of AB is the 
opposite of A'B' (as we have looped once around the origin, the square root is changed into 
its opposite), 

• the contribution of DA is — e'^'^l times the contribution of BC (as we have looped once 
around 1, ln(r — 1) is translated by 2z7r and the orientation of the path is reversed), 

• in the case of v±(s), the contribution of BC is 0(l/|s|) whereas in the case of wo(s) it is 0(1) 
and given by 



VA 



sinh(7r|z|V4) ^/ >p 



|z|V4 



jSinh (7r|zp/4) _ 



r i + z^ 



L- + ^ln|s| 
^ + ^ In 1^1 



.o^(.) = sgn(.)AeW4^,i,,i,(^!L^^+0(^^ 



(28) 

(29) 
(30) 



For later use, we introduce the functions 

a(z)^2Ae--l-lVB -^"M^l^lV4) 

\z\ |z|V4 V 



, h{z) = — e'^^/^v^sinh 



</'(s,^)-y + ^ln|s| 



Using that for a real number u we have r(l + iu) = r(l — iu) and |r(l + iu 
the useful relations 



= -li \ i we get 

|^|2 ^ (2e"|.|V2„2)|a|2, V2|a|2 + |6|2 = V2e"l"l'/4|a|. (31) 



4.2 The third solution 

As an application of Theorem [TT] of the Appendix, let V^{s) and V^(s) denote the two independent 
solutions built in the preceding section through formula and the corresponding A*-^ and A*^ 
contours represented in Figure [2] We set 



\V^{s)\\ \\V^{s)\\ 



(32) 
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and it follows from the preceding result that V'~"^{s) is a solution of (fT5|) . Using ([331), the asymptotics 
of our three (orthonormal!) solutions are 



i^-(^)ii 2|ai \ ^g^^(^), / vn;' 




)+°(r) 



4.3 Scattering matrix 

To summarize, we have 



lim 

s — >-±oo 



%- + ^ln|s| 











e 


4 + ^ln|s| 





( 



where 



1% ^""(-) 



1 / 



± •- 



2 \ay_ 2 |a| +^2 \a\\hy^ ^ 



ie- 

-1 



= M+ 



and the scattering matrix is given by 5 = M+M_ which closes the proof of Theorem [6l 



5 Two-scale Wigner measure and Pseudo Jahn-Teller Hamil- 
tonians 

In this section, we shall prove Theorem [TJ We proceed in two steps: we first prove the propagation 
of two-scale Wigner measure at infinity, then we calculate the transfer's rates at finite distance 
thanks to the scattering result of the preceding section. We shall work in the coordinates (s, 2:, cr, Q 
introduced in Section [3] and then translate the result in the original coordinates {q,t,p,T). For this 
reason, we consider the family (d^) of Theorem [2l It satisfies microlocally in fl, 

op, {Q)v' = ef, 

where for any x G C(^{fl), op, (x)/^ is uniformly bounded in L^(R"^) and where Q is given in (|16p . 
For £ G {0, +1, —1}, we denote by the eigenprojectors associated with the eigenvalues 

of Q. The choice of the labeUing is made so that U J^^^ C {A"*" = 0}. On /, one has zi = 0, 
therefore if (ei, 62, 63) is the canonical basis of C^, 

flO ^ fl- / 62 «) 62 if s > 0, 1 fT+ / 61 (8) 61 if s > 0, , , 

II = 63 (g) 63, II = < ^ „ and = < ^ -r n (33) 

\ 61 (g) 61 if s < 0, \ 62 (g) 62 if s < 0. ^ ^ 
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Let us now investigate the implications of the change of coordinates on the two-scale Wigner mea- 
sures. 



Let V be the two-scale Wigner measure associated with the concentration at the scale ^ of (■0^) on 
/ = {<zAp = 0} with the choice of coordinates on Nil) of Theorem[T] We recall that these coordinates 
(q, r, T]) are such that for p — {q, t,p, r) g /, the class on Np{I) of the vector Sp — {Sq, 5t, Sp, 6t) 
is characterized by the function 

r] ~ Sq A p. 

Let us denote by N{k) the change of coordinates induced on N{I) by k. By the geometric properties 
of two-scale Wigner measures and Egorov Theorem (see [T3i)i the measure P ~ i/ o N{k) is the two 
scale Wigner measure of (w*^) for / affected with the equation 

g{s,z,a,()zi=0 where g o K^'^{q,t,p,T) ^ — — -. (34) 

A(p,t)/x(p,t) 

Besides, D decomposes as 

/> = £>° n° + n+ + i)- n- (35) 

and we are concerned with the link of the traces D^^^. of on S for s > with the traces i'f^ of 
on 5 for s < 0. 

Finally, let us point out that it is not clear that the system ([5]) is well posed in s (indeed, the function 
a depends on the variable a). We address this point in the first subsection, then we successively 
analyze the two scale Wigner measure above |?7| < oo and = oo in the two following subsections. 

5.1 An hyperbolic estimate 

Proposition 8 Let (v^) be a bounded family o/L^(Rs x R^,C'^) satisfying then if (j> ^ C^CR.) 
with (j) = 1 close to 0, there exists So > and £o > such that the family (^(f^')^'^)^^^^, bounded 
m L'^{ILs,L^{Rl)). 

This proposition is proved in the same manner as Proposition 7 in [13j and implies as in [13j that 
one can evaluate the action of pseudodifferential operators on functions (w^) as follows. 

Lemma 9 There exists a neighborhood Qi d D, of such that for all x G C^{Qi), the family 
=opg (x)w^ satisfies: 

fore small enough and for all compact JC C R^, there exists N £ N such that for all a G C(^{IC, C'^'^), 



(op, {a)]f \ v') <C sup sup I a^Sf a(s, z, C) I ds. 



We refer to [13] for the proofs. For the proof of the Proposition, one just need to replace the matrix L 
of [13] by the new definition (fT5|) and for the proof of the Lemma, one has to exchange the role of z 
and C- 

Remark 10 Note that if the symbol of Lemma\^ is of the form a(s,z,cr, C) = ao(z, cr, C)0 (^) 
with r > 0, g (R) and ao G Cjf (R^), then 



(op, {a)v' \lf)\<C 



r-+oo 



ds — > 0. 



Therefore, if p, is a Wigner measure of {v^), we have p ({s = 0}) = 0. 
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5.2 Analysis at finite distance 

We aim at proving Theorem [1] in {\ri\ < +00}, i.e. in {\ri\ < R} for any R > 0. Let us fix i? > 0, 
we use the two- microfocal normal form of Section [221 we localize rj in the ball B = {\r]\ < R} of R 
applying a cutoff function, <j){^) for R > and G C^(R). Then, applying Proposition O we set 
—op/ (Id + ^Ci)v^, and we have 

Va e C;5"(R6 X B), op/ (a) op, {Qo)^ = eF'is, z) (36) 

where Qo is defined in (fT7)) and (F^) is uniformly bounded in L^(Rg2)- 

Obviously, the two-scale Wigner measures of (V^) and (v^) are the same in B and we aim to describe 
the traces on s = 0~ in terms of the traces on s = 0+ of the two-scale semiclassical measures i/ 
of (V^). The function (F^) does not contribute to this description. Indeed, denote by S/s, s') the 
evolution operator associated with the free system (|37p and let (U^) be a solution of this system 
defined more precisely by, 

op, (Qo)C/^ = 0, C/|U-^Lo: (37) 
then we have ^ 

V^s) = U^s) + i [ S/0, t)F''{t)dt. 
Jq 

Since the family {F'^) is uniformly bounded in L^(R^^), we deduce from Holder inequality, 

V^s) = U^s) + 0{VY^\) in L^iRi). 

Therefore, the traces of the two-scale Wigner measure of (C/^) and (V^*^) on s = 0^ are the same. 
In the following, we focus on system ([37)) . which allows to calculate the transfer coefficient in the 
variables (s, z, cr, C) and then, in the variables (g, t,p, r). 

Let us denote as before by i> the two-scale Wigner measure of the family (v^) (and thus of {U^) 
on I X B) for the hypersurface / = {g{s, z, a, ()zi = 0} where g is defined in ([M)) . Because of ([55)1 
and (|35p , one can identify the measures j/" , i/+ and j>~ with the two-scale Wigner measures of the 
components of U^. Set = (J7f , U2, U§), then i/° is the two-scale Wigner measure of (f/f), t'^ of 
(Ul) if s > and of (t/f) if s < and finally, i>~ the measure of (C/l) if s > and of (C/f ) if s < 0. 

The crucial observation is that if one performs a symplectic change of variables (z, i—f (z, C) such 
that zi = a(0, ()zi and if one sets 

U'{s,z) ^ KU'{VIs,V^z) 

where if is a (scalar) Fourier integral operator associated with the symplectic change of coordinates, 
then If^ = {V+,V-,Vo) satisfies the system (fT8| with z — z\. Then, the measure in {s > 0} is 
the two scale Wigner measure of (Kf). By the scattering Theorem [6] and the relation between ujj^ 
in one side and ao, ck+ and a_ on the other side, we obtain 

= (1 - e-^''V2)2 ^+ ^ ^-^ ^ _ ^-.,V2)e-')V2 ,0^ 

on the condition that v^^, and are two by two mutually singular and where the variable 77 
characterize the class of the vector 5p in Np{I) according to 

ry — 5zi. 
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Therefore, 



77 = a(0,C)- 



g(s,2;,cr, C)' 

In view of ([M)) . |A| = a{a, Q ~ ii{<7, ()~^ and of the fact that cr = on /, we obtain 

Studying similarly and i^out^ obtain Theorem [T] in |?7| < +00. 
5.3 Analysis at infinity 

It remains to prove Theorem[l]in {\r]\ = 00}, i.e. that 

^Lt'^H=oo = i>li|^|=oo e {0, +1,-1}. 

Let us argue for the + mode. We follow the strategy of [T3] in Section 7.2 and detail the arguments 
because the fact that we have three modes implies a few complications. Remark [10] shows that 
v+{{s = 0}) = 0. We consider r > 0, oo = ao(z,cr,C) e C^(R^), </> € C^(R), = 1 in a 
neighborhood of 0, (5 < 5o (where (5o is given by Proposition [8]) , i? > 0. We suppose r small enough 
and supp(ao) conveniently chosen so that o(-, rf) is compactly supported in fii for all 77 € R, where Vti 
is the open set of Lemma [HI Then, we have 

(ao , i>+ f l|^|=oo - l|^|=oo) = lim (ao0 (^) , i9.-i>+l|,,|=+oo) • 



We set 



a(s, fj, z, C, ?/) = ao(2, cr, C)0 (^) (y ) (1 " <?^) ( 



The function a is a symbol of A and on the support of a(s, cr, z, C, zi/y^), one has |zi| > so 
that the function an+ is smooth, thus the operator op^^ (afi^) is well defined. Besides, 

( ( - ) , 9si>+l|„|^+oo ) = lim lim lim (op/ (5san+)u'' , 

Therefore, we set 

L, = (op/ (9,an+)w^ , v'^ 
and we shall study as first e ^ 0, then (5 — > 0, i? ^ +00 and finally r 0. 

We remark that 

op/ (a,an+) = - [-a, , op/ (an+)l - op/ (a9,n+). 

Therefore, using the equation, we can write = + + with 

= -i (op/ (a^+)^;^ + i (op/ (a^+)/^ 

- \ 



oPs (ypjT(s,a(cr, O^i)) , op/ (an+] 
(op/ (a9,^+)^;^^;^). 
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We now prove successively that these three terms goes to in the prescribed hmit. We wiU use 
symbohc calculus and for that, we introduce some notations. If q = q{s, z, cr, (, rf), we set 

qe{s, z, CT, C) = 9 ^s, z, cr, C, and ql{s, z, cr, C) = 9 z, ea, eC, 

so that we have 

oVe^ (q) = op, (qe) = opi (ql). 
We observe that 11"'" is homogeneous in the variables (s, a{(7, C)zi), thus for q = an+, we have 

hi 



q,{s,z,a,c) l< (;^) 



where we have used that 



2i,(T,C) :=x/.s2 + a(a,C)2z2>i?Vi (38) 



on the support of q^ and where we implicitely assume < 1. In terms of Weyl-Hormander metric 
(ses sections 18.4, 18.5, 18.6 of [H]), one has g| G S'(l,(7e), where is the metric 

g,^dzl+'^^^±^+e\da'+dC). 

Then 

/ ^\ 2 

9e ^ \/£ 



\ R 

so that the gain of the symbolic calculus is of order -j^. 

• Analysis of L\: Since the function is locally bounded in L^, we consider x G C^(Jli) such 
that X = 1 on the support of a (as before the open set i7i is the one of Lemma IH]). Then the symbolic 
calculus Theorem 18.5.4 of [H] gives 

op/ (afl+) = op/ (afl+x) e opf (an+) op, opi (^S , ^ 

for aU iV e N. Then, the continuity Theorem 18.6.3 of [H] imphes as £ ^ 0, 
(op/ (an+)r I = (op/ (an+) op, (x)r I v^)^o{\). 
Lemma |9] and the boundedness in of op, (x)/"^ gives the existence of C > such that 

(op/ (an+) op, (x)/^' I 

One argues similarly with the other term. 

• Analysis of L^: We first observe that {ypjT{s,a{aX)zi))\ & 5'(0|,5e) where 9 is defined 
in ([38]) and that g| G 5((0|)~^, 5^) for all € N because a is supported in a fixed compact set. 





r+oc 










*(;) 


ds — > 


— 00 
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Therefore, using that VpjT{s,a{a,(^)zi) and 11+ commute, we obtain by symbohc calculus 

1 



OPi 



{iVpMs,a{a,C)zi))l) , opi (q|) 



e opi 



2ie 



i^{(an+)« , {VpMs,a{aX)zi))l}\+l op, (sU^ 



.9e 



It is at this very place that the proof differs from the one of [13] because of the presence of the three 
modes which induces the presence of additional terms. We observe that 

-UVpjT{s,a{aX)zi)t , (an+)»|= ({VpjT{s,a{aX)zi) , {ail+)] 

drjaVpjriO, dc^a{a, , Ozi)i\'^ 



[VpjT{s,a{a,Ozi) , (an+)} 
Similarly, 

i{(an+)« , {VpjT{s,aia,Ozi))l} 

= ff(an+) , VpjT{s,a{aX)zi)} 
Therefore = (op/ {B)v^ , v^) + 0{R-'^) with 
B= \VpjT{s,a{aX)zi) , all+} - \an+ , 

1^ J s.z.cr.C ^ 



^ (d^aIl+VpjT{0,dcMa,Ov) 



Vpjt{s, a 



i'yX)zi)} ^ 



+^d^a (fl+ VpjTiO, %a(a. Or,) + Vpjt(0, %a(a, Qv) n+) 
We write B ^ Bi + B2 + B3 with 

Bi = {VpjT, a}U+ - n+{a, Vpjt} and B^ = a (^{Vpjt. n+} - {n+, Vpjt}) • 



The matrix Bi + ,63 is such that by Lemma [9l we have 

r-+oo / J 



(op/ (Bi+i?3K,.^)|<c|_J(^|0'(^) 



ds. 



The term in ,62 is more complicated and we will cut it into two parts. Indeed, using 

Vpjt{s, aia, C)zi) ^ 

we write 

B2 = a6i (^{n+ -n-,n+} - {n+,n+ -ft-}) +a (^(n+ -n-){6i,n+} - {n+,6i}(n^ 
= ae ({ft+, ft-} - {ft-, ft+}) + a ((ft+ - ft-){^^, ft+} + {e, ft+}(ft+ - ft- 
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We shall prove that B2 is off-diagonal where we say that the matrix D is off-diagonal if H^DH^ = 
for any £ £ {0, -1-1, —1}. We will study those terms together with since dgll'^ is also off-diagonal. 



Set B2 = -82,1 + -82,2 with 

B2,2 = ae ({n+,n-}-{n-,ii+}; 

We remark that the fact that = 5ej,& for aU £, (' £ {0, -1-1,-1} yields 

V<'e {o,+i,-i}, v/teC°°(R^), n^{/i,n±}n^ = 0. (39) 

Therefore, applying ((39|) to h = 6, the matrix B2^i is of the form B2 = aD with D off-diagonal. 



Let us now consider 52,2- Using the general fact A{B,C} - {A, B}C = {AB,C} - {A,BC}, we 
observe that for £^ k, £,k e {0, +1, -1} 

{n^n'^} = {(n^^ ^'=} = ^^{^^^'=}-{^^^'^}^^ (4o) 
-{^^^^} = -{^^ (n^^i ^ jj^l _ Ijjfc n^in^^ 

yielding n°{n+,n-}n" = -frO{n-,fl+}n" = O and 

n~{n+,n"}n" = -n"{n+,n+}n" and n+{n+,n"}n+ = -fi+{n",n"}n+, 
n"{n",n+}n" = -n"{n+,n+}n" and n+{n",n+}n+ = -n+{n",n"}n+. 

We obtain fl^i32,2n^ = for all £ e {0, -1-1,-1} and i?2,2 also is off-diagonal. 

• Analysis of L^: Equation p9)) for h — a implies that dgH^ is off-diagonal. We now consider 
terms of the form 

Ll = (op/ iaD)v' I v') 

for some matrix D homogeneous of degree —1 in the variables (s, a{a, C)zi ) and such that WDW — 
for all i € {0,-1-1,-1}. Without loss of generality we can suppose that D — WDW for some 
£, £' E {0, -|-1, —1} with £ 7^ £' . Then there exists a real number q^^/ {cg^i' £ {±1, ±2}) such that 



[D,VpjT]^ce,e'OD. 



We use this relation to write 

aD -- 



aD 



VpjT 



aD 



CtP" 



-aid + VpjT 



which allows to reuse the equation. Again by Weyl-Hormander calculus with the metric g^, we 
obtain (Vpjt)* G SiOl^g^) and 



aD 



£ S 



1 



■,9e 



(again because of 0| > C^/eR on the support of a|). Therefore, we have 



op/ (aD) £ opg {-a + Vpjt) op/ 



aD 



aD 



+ op/ ( ds 



aD 



op^ (-0- 4- Vpjt) 



OPi ( 5 ( 
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whence, using the equation for v'^ 



r I v' 



Arguing as before, one gets as e goes to then R to +00, 



I v' 



Observing that 



8 

cr.C 



aD 



< 



we obtain by Lemma [9] 



1 \ /•+00 
\Ll\<o{l) + 0( — ]+Ce 



(s2 + ei?2)3/2 



,<o(l) + - 



(s2 + ei?2)3/: 



which finishes the proof. 
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A Appendix: Complex wedge product and solutions of ODEs 

We consider (7i,(-|-)) an Hilbert space of dimension 3 over C and A{s) a continuous family of 
selfadjoint endomorphisms on Ti.. We consider the Schrodinger system 



- id,V = A{s)V. 



(41) 



Let [•,•,•] denote a non trivial alternating trilinear form on Ti. and let ui and U2 be two vectors of 
H. We define the wedge product of ui and U2 with respect to [•,•,•] to be the unique vector ui A U2 
(through Riesz representation theorem) such that 



y-w e H, [ui,U2,w] = {ui A U2\w). 
Then, we have the two following points: 
1. If A is an endomorphism of Ti., we have 

[Aui, 1*2,^3] + [UI,AU2,U3] + [mi,U2, Aus] = (TtA)[ui,U2,U3]. 



(42) 



(43) 



2. If u (respectively v) has coordinates (mi, 1*2,^3) (respectively (wi,W2,i'3)) in some orthonormal 
basis (ei, 62, 63) of C"^ then the coordinates oi u A v are 



[ei, 62, e3](w2W3 - U3V2,U3Vl - UlV3,UiV2 - U2V1). 

Theorem 11 Let Vi{s) and V2{s) be two solutions of the system Then 

6 -^=0 Vi{s) A V2{s) 



(44) 



(45) 



is also a solution of In particular, if Vi and V2 are linearly independent, we have a basis of 

solutions. 
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Note that this is in fact a corollary of Liouville theorem on Wronskian combined with the fact that 
the equation is norm preserving. 

Proof: If Vi and V2 are linearly dependent, the wedge product is zero and is certainly a solution. 
Thus, we can assume that Vi and V2 are linearly independent so that (Vi(s), ¥2(3), Vi{s) A V2(s)) is 
an instantaneous basis of H. Removing the phase, it is enough to show that W{s) — Vi{s) A ^2(5) 
is a solution of 

-idsW = [A{s) ~ TTA{s)ld]W 
which we will show to hold in the preceding instantaneous basis. We first have 

{V,{s)\{~zdsW){s)) = {-zdsV,\W) - {AV,\W) = (VMW) = - TrA(s)Id]M^(s)) 

where we have used that Vj and W are orthogonal, the equation and the self-adjointness of A{s). 
Similarly, using and ([33]), we obtain 

{W{s)\{-id,W)is)) = ~{W\{-idsVi)AV2)-{W\ViA{-id,V2)) 

= --{W\AViAV2) ~{W\ViAAV2) 

= -[AVi,V2,W]-[Vi,AV2,W] 

= [V,,V2,AW]-{TtA)[VuV2,W] 

= {[A- TrAId] 14^1^1 A V2) 

= {W{s)\[A{s)-TrA{s)ld]W{s)). 

□ 
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